APS/123-QED 



o 

(N 

1—5 



I 



cd 



I 

O 

o 



> 

00 

rn 

O 



X 



Tunneling between helical edge states through extended contacts 

G. Dolccttoi^2,3^ g Barbarinoi, D. Ferraroi-^'^^ N. Magnoli^^^, M. Sassetti^'^ 

Dipartimento di Fisica, Umversita di Genova, Via Dodecaneso 33, 16146, Genova, Italy. 

^ CNR-SPIN, Via Dodecaneso 33, 16146, Genova, Italy. 

^ INFN, Via Dodecaneso 33, 16146, Genova, Italy. 

(Dated: February 27, 2013) 

We consider a quantum spin Hall system in a two-terminal setup, with an extended tunneling contact 
connecting upper and lower edges. We analyze the effects of this geometry on the backscattering 
current as a function of voltage, temperature, and strength of the electron interactions. We find that 
this configuration may be useful to confirm the helical nature of the edge states and to extract their 
propagation velocity. By comparing with the usual quantum point contact geometry, we observe 
that the power-law behaviors predicted for the backscattering current and the linear conductance 
are recovered for low enough energies, while different power-laws also emerge at higher energies. 

PACS numbers: 73.23.-b, 71.10.Pm, 73.43.-f 



I. INTRODUCTION 

Since the discovery of the Quantum Hall Effect 
(QHE)i, the condensed matter community devoted great 
efforts in finding other topological states of matter in 
which fundamental physical properties are insensitive to 
smooth changes in material parameters and can be mod- 
ified only through quantum phase transitions. In recent 
years, a new class of these peculiar systems have been 
experimentally observed: the topological insulators^i^. 
Their main characteristics are the presence of a gap in 
the bulk, analogous to the one of the ordinary band 
insulators, and gapless edge states protected by time- 
reversal symmetry. In two spatial dimensions they are 
the first realization of the Quantum Spin Hall Effect 
(QSHE), theoretically predicted in graphene with spin- 
orbit intcractioriii^, in strained semiconductors^ and in 
Mercury- TcUuride quantum wells^*^. The edge states of 
the QSHE are helical^, namely their electrons have spin 
direction and momentum locked each other. In presence 
of intra-edge interactions they can be described in terms 
of a helical Luttinger liquid^. The experimental mea- 
surement of non-local transport in multi-terminal setups, 
according with the prediction of the Landauer-Buttiker 
theoryiSj represented an important test of the existence 
of helical edge statesii. 

The fast technical developments in this field will make 
shortly possible to realize interesting experimental ge- 
ometries, like the Quantum Point Contact (QPC)^^"— 
that already revealed extremely useful to extract in- 
formation on the edge properties in the fractional 
QHEi^"— . Various theoretical proposals have investi- 
gated this geometry focusing on both the two-terminal2i 
and four-terminal^^Ti^i setups. Possible interference 
expcrimcnt a^^'^^ , as well as quantum pumps^, involv- 
ing two point contacts have also been considered. 
The possibility offered by the Mercury- Telluride quan- 
tum wells to realize a QPC by means of electrostatic 
gates or, more realistically, by etching the sample in 
the desired shape makes possible to have a great con- 
trol on the geometry and allows to study the evolution 



of the transport properties as a function of the constric- 
tion geometrical parameters. An analysis of the effects of 
extended contacts^^Ti^i on the transport properties have 
been already addressed for the QHE showing deviations 
from the standard power-law behavior of the current as a 
function of the voltage at zero temperature. Finite tem- 
perature effects were also considered for composite frac- 
tional QH systems^, demonstrating that extended con- 
tacts may provide information about the neutral mode 
propagation velocity along the edge, provided that it is 
very small with respect to the one of the charged mode. 
In this paper we propose to investigate the extended con- 
tact geometry for the helical edge states of the QSHE by 
properly taking into account the role played by interac- 
tions. We will evaluate the backscattering current as a 
function of voltage and temperature. We will demon- 
strate that all the deviations with respect to the point- 
like case can be included in a modulating function. We 
will demonstrate that, at low enough temperatures, a 
peak appears in the differential conductance, which pro- 
vides evidence of the helical nature of the edge states and 
gives information about the propagation velocity of the 
edge modes. At low energies the backscattering current 
and the linear conductance arc described by the same 
power-law behaviors predicted for the QPC geometry. 
Even more interestingly, power-laws are recovered also 
at higher energies, but with different exponents. 
The paper is divided as follows. In Sec. |ll]we recall the 
main results of the helical Luttinger liquid description of 
edge states of a QSH system. In Sec. IIIII we analyze the 
extended contact geometry introducing the modulating 
function both in the non-interacting and in the interact- 
ing case. Sec. IIVI contains the main results on transport 
properties. Sec. |V]is devoted to the conclusions. 



II. MODEL 

We consider a QSH insulator with one Kramers dou- 
blet of helical edge states in the two terminal configura- 
tion (see Fig. [T]). On the upper edge (1) one has right- 



moving spin up and left-moving spin down electrons, on 

the lower edge (2) the opposite. 

The corresponding free Hamiltonians aro^ii^ (fi, = 1) 



Hi 



(2) 



-ivp I dx 



^kT(;)'^^'^«^t(;) 



■ V'l,;(t)^^V^L,;(T) 



(1) 
where ipB.,i; ('/'i.t) annihilates right (left)-moving electron 
with spin up, and analogous for the spin down, and vp is 
the Fermi velocity, estimated^i'^ about 5 • 10^ m/s. For 
sake of simplicity we assume infinite edges, even if a more 
realistic description based on finite length edges coupled 
to non-interacting leads can also be considered^^i^. This 
so called g{x) model^^"— reveals crucial in order to re- 
cover the proper quantization of the conductance of one 
dimensional channels and leads to finite length correc- 
tions to physical quantities, that however are not crucial 
in the considered setup^^. 

Concerning interactions, we consider terms which pre- 
serve time-reversal symmetry near the Fermi surface for 
a single Kramers doublet of helical edge states^. They 
are a subset of all possible contributions analyzed by the 
so called g-holog y^'^'i'^^ represented by the dispersive 

Hd = ,g2_L j dx ('/'l{,^'/'fl.T^U'/'i4 + V'l.^.V'L.t^k^'/'fl.i) 

(2) 
and the forward scattering 



Hf 
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— ^ / dxV'l^^Va^o-V'l.^V'a.^T. (3) 



Note that possible Umklapp terms, which are impor- 
tant only at certain commensurate fillings^, are here ne- 
glected. 

The bosonized procedure of the Luttinger liquid allows 
to write the electronic field operator in the form^ 



i^R/bA^) 



T. 



R/L,cr ^±ikFX^- 



■iV2TTipii/L^„{x) 



(4) 



with ipR/L.aix) a bosonic field (<t =t,i), ■FR/L,a the 
Klein factor, necessary to give the proper commutation 
relation between electrons belonging to different edges, a 
a finite length cut-off and kp the Fermi momentum. The 
bosonic field ^PR/L.a{x) is related to the electron density 
through PR/L,aix) = T^d^,ipii/L^cr{x). According to 
the standard bosonization procedure^^i^ the interaction 
terms in Eqs. ([2])- ([3]) are quadratic in the electron den- 
sity. 

Introducing the helical edge basis on the upper and lower 
edge^^ 

¥'1(2) (a;) = -y= [ipL,Ui)i^) - fR,W)i^)] ' (5) 

with their canonical conjugates 

^1(2) (a;) = -y| [<PL,ui)i^) + "PR.imi^)] ' (6) 



the total Hamiltonian H — Hi + H2 + Hd + Hf can be 
recast in the bosonized form^iiS^ 



H ^ 



-y 
2 ^ 

i=l,2 



dx 



-^ {dxfif 



K {d,0. 



(7) 



Here, K ~ \/ nZl'l 1 f'^" 1 ^J^ is the interaction parameter 



2iT'UF+g4|i+g2 



and V 



Vp 



1 
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the rcnormalized 



2ttvf j \ 2ttvp 

velocity. For Coulomb repulsion g^ = g2± and therefore 
V — vp / K. In the following we will assume this condi- 
tion, despite other possible interactions can be straight- 
forwardly taken into account. 



III. EXTENDED CONTACT 

In presence of an external voltage V , right (left) 
-moving electrons feel a chemical potential /i^ (/^a), 
with /ii — /ifl = eV . Spatial separation prevents 
electron tunneling between edges leading to conductance 
quantizations Q = S- in order to study tunneling effects 
the system is pinched by means of a gate voltage^i or, 
more realistically, by etching the sample^ creating a 
tunneling regioni^. 

Previous theoretical works have studied this 
configurationi^i^"— i2^, both in two-terminal and in 
four-terminal setups, assuming a point-like tunneling. In 
what follows, we will generalize this assumption, taking 
into account the possibility of tunneling events occurring 
in an extended region (see Fig. [ij. Our aim is to 
investigate the effects induced by a long contact on the 
backscattering current. The backscattering Hamiltonian 




Figure 1. (Color online) Extended contact geometry for a 
quantum spin Hall system with one Kramers doublet of heli- 
cal edge states. The full (dashed) lines represent hehcal edge 
states carrying electrons with spin up (down). Right (left)- 
moving electrons are in equilibrium with the left (right) con- 
tact at chemical potential fiL (ma)- The black arrow rep- 
resents a possible spin-conserving electron tunneling event 
through the extended region. 



connecting the two helical edge states is represented by 



Hv 



dx dy 



Yl ^^,v^R.A^)^LAy) 



=t.l 



+ h.c., (8) 



with A^_y the tunneling amplitude in which a left-moving 
electron is destroyed in y on one edge and recreated as a 
right-moving electron in x on the other one. A reasonable 
choice for A^.y is to assume a separable forro^S 



Ax.y^Aofi{\x + y\)fci\x~y\). 



(9) 



The function /;, indicated as lateral contribution, spec- 
ifies the average location of the tunneling events^^""— 
while fc, dubbed crossed, allows to take into account 
non perfectly vertical events^S. This assumption is rea- 
sonable for smooth tunneling junctions. Both functions 
are maximal around zero and decrease by increasing 
their arguments. With this requirement, the longer 
the tunneling path the smaller the corresponding local 
amplitude. 

Note that Eq. ^ describes spin-conserving tun- 
neling processes only, since spin-flipping tunneling 
terms give no contribution in our two-terminal 
setup ^^'^^ . Furthermore we neglect tunneling of 



either charged 
ful (- c 



(^ 



cos 



2) ) 



^ (<y3i + <^2)J ) or spin- 
particle pairs, although 



7^ (Si --2) 

for strong enough electron interactions they could 
compete with single-particle tunneling processes 



(^ 



ib^i'Pi + 'p^] 



2V^ 



1 

2V^ 



^(^1-^2) )^^. Note 



that all these processes are irrelevantii^, in the RG sense, 
for 0.5 < K < 2. We limit our analysis to repulsive 
interaction 0.5 < K < 1, and we treat the tunneling 
current as a small perturbation. 

The tunneling Hamiltonian in Eq. (|5]) induces no net 
charge transfer between the two edges, but leads to a net 
spin tunneling current. The corresponding spin current 
operator is 

k = -^ Yl / ^^ ^2/ ^x,y1pRAx)-)pL,a{y) + h.C, (10) 

according to the requirement of absence of spin flipping 
and multiple-particles contributions. In linear response 
approximation in the tunneling hamiltonian, the station- 
ary expectation value of the spin current in Eq. (|10|) can 
be written in terms of the tunneling rates Ti^^^^ji^a and 

A R,a—fL,a SS 



(/s 



4i: 



■ L,<T^rR,a 



R,a^>L,a\ 



(11) 



a=t,i 



Note that the functional dependence of rates and other 
physical quantities from bias and temperature is under- 
stood for notational convenience. 

One can easily realize that this spin tunneling current 
is responsible for a reduction of the net current flowing 



from one lead to the other— i2^, i.e. (/) = ^V — (/bs), 
with (/bs) the backscattering current, related to (Ig) by 



/bs> = 2e(/s 



(12) 



We can thus measure the spin tunneling current by mea- 
suring the ordinary backscattering curren1j2i. 
By taking into account the spin independence of the tun- 
neling rates and by considering the detailed balance re- 
lation rfi,^_,L,a = e-f^^^TL^a^R^a, {P = l/ksT thc in- 
verse temperature) one has 



(/] 



Bs; 



2e(l 



^-PeV\ 



■ L,t^i?,t- 



(13) 



According to Eq. p^ . we can consider only the tunnel- 
ing rate T = Tl,-\^r„^ given by 



r = dx dy dx' dy' Ax^yA'^ry 

dte'''''G>{y'-y,t)G<{x'~x,t), 



with 



^R/bi^^^) 



'R/L 



„WR/L{x,t) 



'-'r/l(^^*) 



2'Ka 



2'Ka 



(14) 

(15) 
(16) 



the greater and lesser electron Green's functions associ- 
ated to the right {R) and left (L) movers. The corre- 
sponding bosonic Green's functions are 

WR/L{x,t) = 2TT{<fn/L,a{x,t)ipR/L,a{0,0)) 

-2^((^flyi^,(0,0)(^flyi,,(0,0)). (17) 

They do not depend on spin and can be written in terms 
of the chiral ones 'W±{x, t) 

Wfl(a^,0 = c^^W+(x,t) + c(^%_(.T,i) (18) 
WL{x,t) = c^i,^W+{x,t) + c^+^W-{x,t), (19) 



with 



and 



W±(a;,i)=W iT- (20) 



W{t) = In 



{}+Tk-^i) 



Y2 



(l + ^)(l + ^'^ct) 



(21) 



.(±) 



Here, T{x) is the Euler Gamma function, c\^ 

■J ( yK ± -^ 1 are the interaction dependent tunneling 

coefficients and ujc = v/a thc energy bandwidth. By re- 
placing the above expressions into Eq. (J14p one obtains 



A A*, , 



(27ra)2 



"^7 



Tk= \ dx dy dx' dy' /' ^f ^k.i^y'^yW-^ / ^^ ^.eVt^c'~-'wit-^^)+c^-'wit+^-^)+c^-'w(t-y-^)+c<^'wit+y-^) ^ 



(22) 



where we explicitly indicate the dependence on the inter- 
action parameter K . 



I 

This can be done regardless of the form of the tunneling 
amplitude but, as we will see, the separability assumption 
In what follows we will first analyze the non-interacting of Eq. ([9]) allows to give a closed form for the modulating 
case, which can be thought as a superposition of two function. From Eq. ([T51) and Eq. ([25]) follows that 
independent integer QH systems subjected to opposite 
magnetic fields. Later we will address the case of inter- 
acting helical edge states. 



(/Bs) = Aix(4r"*)). 



A. Non-interacting helical edge states 

In the non-interacting case {K ~ 1), one has c)^Li = 1 
and Cj^Li = 0, and Eq. (P^ reduces to 

' dt e-^*e^(*-'^)+^(*+'^^ (23) 






where we introduced the short hand notation dx = dx ■ 
dx' , dy = dy ■ dy' . In terms of the new variables^ r = 



Tl 



ds dy ^^;y^*-''y' e'M^'-^)+'-iy'"y)] 

^ (27ra)2 



X dr e 



ieVr [W(r-^)+W{T+^)] 



(24) 



with k± ~ kp ± eV/2v. This can be further expressed as 



Fi = / dx dy 



(27ra)2 
where 

Fg{z,uj)^ JdTe'"^Pg(^T 

and Pg{t) = eff^(*) (cf. Eq. ([21 



yL^^[k+{x'~x)+k^{y'-y)]p^^^^^y^^ 



Pgir 



(25) 
(26) 



The separability assumption in Eq. (jH]) allows to factorizc 
the tunneling amplitude as 



I lAor 

(27ra)2 



dy cos 



eV 



iy' - y) 



U\2y\)UW\) 



X / dxcos[2kF {x' - x)] fi{\2x\)fi{\2x'\)Fi(x' - x,eV). 

(27) 

To better characterize the effects of the extended contact 
geometry it is useful to represent Fi in terms of the point 
contact rate r^^°™ ' as 



Fi = Ai X F^^'""*^ 



(28) 



'Bs /• (29) 

For any interaction K, the point-like current is given by^ 

|2 



,j(point)^ __ ,^^^ _ ^^peV^ 1^0 



0=2e(l-e-'^^^) 



.„ ,^P2dAeV) (30) 

(27ra)^ 



with dK = Cx 



„(-) 



= i (K + i) so that dx = 1 in 



K 



the non-interacting case. The function 



Pg{u;) = dt e^-'Pgit) 



(31) 



has the following forrr^ for energies lower than the band- 
width OJn 



.9-1 



P^iE) 



fe(£) ' (^(E) (T^O) 



r(s) 

2tt ] e~a~ n 



2 ' 2ir ' 2 "■" ' 27r 



(T^O) 

(32) 
with d{x) the Heaviside step function and B [x, y\ the 
Euler Beta function. 

The modulating function Ai in Eq. (|28p represents the 
influence of the extended region and is given by 



Xi = 4 dy cos 



eV 



(y' - y) 



U\2y\)U\2y'\) 



X dx cos[2kF {x' - x)] fi{\2x\)fi{\2x'\) 
Fi(x' -x,eV) 

X ^t: -. 

P2{eV) 



(33) 



It can be written as a product of crossed and lateral con- 
tribution Ai = XiX\; with 



X'{^2 i dy cos 



eV 



{y - y) 



U\2y\)U\2y'\) (34) 



A'l = 2 / dx cos [2kF [x' - x)] /,(|2.t|)/z(|2x'|) 
FAx' ~x,eV) 

X ^t: -. 

P2{eV) 



(35) 



Notice that, while A^ depends on the crossed contribution 
fc only, X{ contains also the electronic Green's functions 
through Fi. 



In order to perform an analysis of the extended contact, 
we consider a separable gaussian form^Si^ 



A 



An 



x,y 



27tUi 



(36) 



The parameter ^; is related to the extension of the con- 
tact, while ^c allows to take into account non perfectly 
vertical events. In this sense a realistic assumption for 
modeling an extended contact is ^c "^ 6- Note that in 
the limits £,c,i -^ Owe recover the point- like tunneling am- 
plitude A^.y -^ AoS{x)d{y), so that (Ibs) ^ (-^55™*^)- 
By replacing the gaussian expression into Eqs. (|34p - ([55)) 
one obtains 



X{ = 



1 f ic<!v \^ 
e 2\ ^ ) 



(37) 
dxe-"^ cos {2kF^ix) £lte^l£Zl. (38) 



/2-K J •. ' ^- , P2ieV) 

By exploiting the convolution properties 



Fg{z,Uj) 



1 

2^ 



dE e 



i^E 



Pn 



E\P„ 



E 



.2 J " \2 

the tunneling amplitude can be written in the form 

Ai = 



(39) 



^-K«c^)=-2(^.6)^ j ^e~'(^'^r cosh (4fc^e?f ) 



P-zieV) 



(40) 
I 



This result is valid also at finite temperature and ex- 
tends what done in Ref. io| for the QHE at T = 0. Note 
that the crossed contribution to the modulating func- 
tion comes into play only at high bias voltage. For an 
extended contact with length ^ (0.1 -^ l)/-im, one has 
^i ^ (O.l-^l)^m and (c < 6, e.g. Ic ^ 10 nm. With this 
assumption the crossed contribution is crucial only for 
relatively high bias > 0.1 V, not considered here. This 
fact allows to choose A^ w 1 and to focus only on the lat- 
eral contribution which, as we will see in the following, 
shows strong modifications with respect to the point-like 
case also at low bias. 



B. Interacting helical edge states 



Starting from the general expression in Eq. (p2|) and 
proceeding as in the previous section, one can express the 
interacting modulating function as {K ^ 1) 



\K 



_ f dEidE2dEz -^[^{eV-2E2-2E'i)Y-i^^^(eV-2Ei-2E2-2kFv)^ 



,PA+AEi)PA-^{E2)PA-AE^)PA+){eV- E E,) 



P2dAeV) 



4=1,2,3 



(41) 



Due to the natural constraints imposed by the functional 
form of P{E) in Eq. (|32p it is possible to neglect the 
crossed contribution, present in the first gaussian term, 
as far as eV, ksT <^ v/ic- Under this condition and 
noting that 



dp 

— Pg, iE)Pg, iu;-E)= Pg,+g, (w), (42) 

ZTT 



Eq. (|4T|) becomes 



^ ^ „-2af I "-^ ^-2{Kai- 



dE _2 



cosh ( AKa^ — j 



Pa^{^ + E)Pa,{^^E) 
P2dAeV) 



(43) 



Here, we introduced the Fermi energy ep = kpvp and 
the dimensionless parameter ai = kp^i- The modulating 
function thus depends on the length of the contact ^; and 
on the Fermi momentum only through their product. By 



inserting Eq. (|5^ in Eq. (|15)) one has 

,„, . ffig£_!l /,„-*'-' ^->»h (2Ka;tsl. 

87r^r^(rfif) J \ ep 

X 6[7+,+ (x),7+,-(x)]6[7-,+ (.T),7-,-(2:)] (44) 



with (77, 7/ = ±) 



, . dx i ( eV , 



(45) 



To conclude we observe that also in the interacting case 
the backscattering current can be written as 

(/bs (V^, T)) = \k {V, T) X (47"*) {V, T)) (46) 

with (/b^°'" \V,T)) given in Eq. ([29]) and where we ex- 
plicitly reintroduced the dependence on bias and temper- 
ature. Note that for ai = 0, Eq. ([^^ reduces to A^- = 1, 
and the point-like tunneling case is recovered. 



IV. RESULTS 

Since the modulating function depends on bias and 
temperature, it will influence the behavior of transport 
properties with respect to the point-like tunneling case. 
It is then useful to investigate it in details. Fig. [2] shows 
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Figure 2. (Color online) Modulating function as a func- 
tion of (a) bias V (in units of (.f /e) at low temperature 
{ksT — 10~^eF) and (b) temperature T (in units of ep /ks) 
at low bias {eV = 10~ ef), for different lengths of the con- 
tact: a; = I (long dashed red), 2 (dashed green), 5 (short 
dashed blue). Note that the behavior at low temperature in 
(a) is indistinguishable from the T = case. This comment 
holds as well for panel (b) between low V and V = Q. Other 
parameters: K = 0.75. 



Xk as a function of voltages (a) or temperatures (b). 
Fig. [2l^a) presents a maximum sA V ~ V = 2eF/eA", 
becoming more and more pronounced by increasing a;, 
that is the length of the contact. In the limit a/ — )- 
it is washed out and \k[V,T) -^ I. As already noted 
for QHIiSSi, this maximum is determined by the two 
phases that control tunneling, one set by the Fermi 
momentum {2kpx) and the other by the voltage drop 
(eVt). The peak occurs when the two phases are equal: 
eV = 2kFx/t = 2kFV = ItFJK. 

A maximum is present also in Fig. El^b), but it originates 
from a dephasing mechanism, induced by finite temper- 
ature, similar to what was found in interferometric ge- 
ometries with two or several QPCs, both in QH^ and 
in QSH systems^^, where the dephasing was depending 
on the distance among the QPCs. The extended con- 
tact geometry can be seen indeed as an infinite series of 
QPCs with different tunneling amplitudes, with infinites- 
imal distance dx between them, and the backscattering 
current is now given by integrating over the contact re- 
gion. For all interaction strengths 0.5 < X < 1 we find 
the maximum at a position T of the order of ep jkQ^ van- 
ishing as a/ — )■ 0, reproducing in this case the point-like 
regime with Aif(y, T) — )• 1. 

Note that for vanishing bias and temperature the modu- 
lating function is exponentially suppressed by the length 
of the contact, namely Ak(^ = 0,T = 0) = g-^"?. 
We can also study the asymptotic behavior of \k at low 
bias or low temperatures. Introducing the energy scales 



eVai = eF/{Kai) and ksTai = €F/{Kai) one finds 

constant V ^ Vq, 



XK{V,T^eV/kB 
and 



V- 



v-v:^ K., 



(47) 



\ (M ^^ 1 rri rpN J constant T -^T^i ,.„^ 

\K\y <-kBTIe,T)^ i. y_i T-f^T ' ^ ' 

Fig. [3] shows the differential conductance G{V, T) = 
d{lBs{V,T)) / dV as a function of bias (a) and the lin- 
ear conductance G{T) = G{V = 0, T) as a function of 
temperature (b). They both show a peaked structure, in 
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Figure 3. (Color online) (a) Differential conductance as a 
function of bias V (in units of ef/g) at low temperature 
{ksT — IO^^ef) and (b) linear conductance as a func- 
tion of the temperature T (in units of ep/kB), for different 
lengths of the contact: ai — 1 (long dashed red), 2 (dashed 
green), 5 (short dashed blue). Units of the conductance: 
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(27ra)- 
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Other parameters: K = 0.75. 



contrast to the point-like case, reminiscent of the form of 
Xk (see Fig. ^. 

More quantitatively, focusing on a given length, we can 
study the dependence on interactions. Fig. |4] shows the 
differential conductance as a function of bias, varying the 
electron interaction. The conductance shows a peak at 
V ^ V, which depends on the velocity of the excitations 
{V = 2kFv/e). Thanks to this behavior, we argue that 
an extended contact geometry could be fruitful to extract 
information about the velocity of the excitation modes 
along the edges, by experimentally measuring the peak 
of the conductance, varying the Fermi energy and the 
bias voltage^iii. Furthermore, it must be stressed that 
in presence of an ordinary Luttinger liquid we should 
expect two different peaks, as a consequence of the spin- 
charge separation, which leads to two different propaga- 
tion velocities, one for the charge modes and one for the 
spin modes^ii^iilii^. The single-peak structure of Fig. 
m instead, provides evidence of the close connection be- 
tween spin and charge typical of the helical edge states 
of QSHE, where these degrees of freedom are locked each 
other and propagate with the same velocity. 
We remark that, as expected, the peak in the differen- 
tial conductance is reduced by increasing temperature 
and finally washed out for temperatures 2'nkBT ~ eV , 
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Figure 4. (Color online) Differential conductance as a function 
of bias V (in units of ep/e) for different interaction strengths; 
K = 1 (long dashed red), 0.75 (dashed green), 0.5 (short 
dashed blue). Note that the conductance is plotted in unity 
of Go as in Fig. [S] which depends on K and thus not allow for 
a direct comparison on the size between the different curves. 
Other parameters: ai = 5; ksT — 10~^eF. 




Figure 5. (Color online) Differential conductance as a function 
of bias V (in units of ef/e) for different temperatures (in units 
of eF/fcs): T = 0.1 (long dashed red), T = 0.5 (dashed green) , 
T = 2 (short dashed blue). Units of Go as in Fig. El Other 
parameters: ai =5; K = 0.75. 



as shown in Fig. [5] 

Information about velocity are not the only ones that 
can be extracted by means of this setup. Theoretical 
works concerning point-like tunneling predict power-law 
behaviors for curren t^^'^^ 



{lil°'''*\v,T«eV/kB))^V 



2rfK-l 



/^(point) f^\ Y^' 



2d,. 



(49) 



(50) 



Despite these trends arc here no longer valid, they still 
survive at bias or temperatures lower enough, namely for 
V <C Vai or T ^ Tq, respectively, as shown in Fig. [51 
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Figure 6. (Color online) Log-Log plot (a) of the backscat- 
tering current (in units of Jo = (eF/e)Go) as a function 
of the bias voltage V (in units of tp/e) at low tempera- 
ture (fcsT = 10~^eF) (long dashed red curve) and (b) of 
the linear conductance (in units of Go) as a function of 
temperature (in units of ep/kB) (long dashed red curve). 
Other parameters: ai = 1, K = 0.75. Straight lines rep- 
resent the asymptotic power-law behavior with exponent (a) 
2dK - 1 = 13/12 (dashed green line) and 2dK - 2 = 1/12 
(short dashed blue line) and (b) 2dK — 2 = 1/12 (dashed 
green line) and 2dK — 3 — —11/12 (short dashed blue line). 



Interestingly, by increasing energies, new power-law be- 
haviors arc recovered, however, with different exponents 



and 



(V^-F»F„,) 

(51) 



G{T) - T 



■2(iK-3 



{T-f:$>Ta„). 



(52) 



This is a consequence of the asymptotic behavior of the 
modulating function (cf. Eqs. (|T7)) - (|i5)) ). 
It is worth noting that the effective visibility of these high 
energy power-laws crucially depends on the Fermi energy 
ep of the system, that can be easily tuned experimentally 
by means of an external gate^, and the natural cut-off 
energy Wc of the theory. The latter can be reasonably 
identified as the energy at which additional bulk effects 
have to be taken into account, thus the presented helical 
Luttingcr liquid picture holds for energies lower than Uc. 



V. CONCLUSIONS 

We proposed a model for an extended tunneling 
through contact region in QSH system. We demonstrated 
that it is possible to take into account the extended na- 
ture of the contact through a modulating function, which 
rcnormalizcs the transport properties of the point-like 
case. 

We showed that, due to the extended nature of the con- 
tact and for low enough temperatures, the differential 
conductance shows a pronounced peak that can be used 
to extract information about the propagation velocity of 
the excitations along the edge. The presence of a unique 
peak is a signature of the helical nature of the edge states 
in QSHE. 



We analyzed the backscattcring current in the low tem- 
perature regime and the linear conductance, showing that 
the power-law behaviors predicted in the point-like case 
survive at progressively lower energies by increasing the 
length of the contact. Remarkably enough, new power- 
laws emerge also at higher energies, but with different 
exponents. 
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